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This contribution describes an approach to the problem of resolution and resolvability in scattering 
methods (e.g. X-ray diffraction, electron microscopy) in the presence of series termination and 
blurring. One-, two- and three-dimensional cases are considered separately. Formulas relating the 
effects of nominal resolution and blurring to peak resolvability are derived and analysed. We show 
that both blurring and series termination widen point source peaks thus reducing their resolvability. 
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INTRODUCTION 


Macromolecular crystallography (MX) and 
single particle cryo electron microscopy (cryo- 
EM) are two widely used scattering methods used 
to derive atomic models of biological molecules 
(Rupp, 2010; Frank, 2005). Both techniques pro- 
duce density maps that are interpreted as atomic 
models (Emsley and Cowtan, 2006). Derived ato- 
mic models are deposited to the Protein Data Bank 
(PDB) via which they are made available to the 
community free of charge (Berman et al., 2002). 
Maps can be represented equivalently either in ima- 
ge space as a function sampled on a regular three- 
dimensional grid points or as Fourier coefficients of 
these maps. In crystallography observations are re- 
lated to the Fourier coefficients (Harker and Kasper, 
1948), whereas in cryo-EM usually observations are 
made for images themselves (Frank, 2005). Details 
that can be seen in these maps and therefore accu- 
racy of atomic models depend on the amount and 
quality Fourier coefficients. 

There are some confusion and controversy on 
the definition and use of resolution concept: cur- 
rent definition only uses the highest frequency 
terms of the Fourier coefficients used for map cal- 
culations (Wlodawer et al., 2017). This definition 
does not account for the noise in the data, absolute 


and relative mobility of the objects under study. It 
only states the radius of a sphere within which ob- 
servations reside. Although it does seem that in 
crystallography, since observations are directly re- 
lated to the Fourier coefficients, the highest obser- 
ved frequency term should define the resolution 
provided that all Fourier coefficients have been 
measured with high accuracy. This definition has 
been successfully used for last 100 years. However 
now increasingly more and more noisy data are 
used for structure determination and old definition 
is no longer applicable: this definition only defines 
the radius in the Fourier space where data reside. 
Another problem arises when the resolution of 
structure used as an indication of their accuracies 
and different structures are compared with referen- 
ce to resolution; e.g. for selection of the “best” ato- 
mic models, for comparison of accuracy of atomic 
models to typical models in the PDB correspon- 
ding to the same resolution. 

In cryo-EM resolution is defined by analogy 
to that in crystallography: the highest frequency of 
observed Fourier coefficient are defined for which 
Fourier shell correlation (fsc) is more than 0.143 
where fsc is calculated in narrow shells between 
Fourier coefficients of two independently reconst- 
ructed maps (Rosenthal and Henderson, 2003; 
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Scheres, 2012). Both of these definitions, to a cer- 
tain degree, make use of the signal to noise ratio. 
Information about the signal-noise ratio should al- 
so be used in structure determination and refine- 
ment (Murshudov, 2016). Usage of the highest ob- 
served significant frequency terms as a definition 
of the resolution does not account for the behaviour 
of the noise and the signal over the Fourier space: 
it is time to make these definitions more precise 
that accounts for the quality of all available infor- 
mation in the data. 

Moreover, both of these definitions do not 
fully account for the signal the noise level in the 
data, they also do not use the fact that molecules 
oscillate resulting in the blurred image and there- 
fore reducing the visibility of them in the map af- 
fecting different parts of the density differently in 
accordance with varying mobility of molecules. 

In this paper, we will address this problem 
theoretically with some certain simplifying as- 
sumptions that make calculations manageable. He- 
re we consider one-, two- and three-dimensional 
cases and show that objective measure can be cal- 
culated to measure resolvability that would reflect 
the information contained in the data about the ob- 
ject under study. 


1. Resolution and Resolvability 

Let us assume that we have observations in the 
Fourier space up to a nominal resolution Snom. Our 
observations are Fourier coefficients that are sum 
of two components Fon5(S)=Ftrue(S) + Froise(s), 
where Fops(S) is observed, Firue(S) is the true” 
and Fy oise(S) is the noise Fourier coefficients. We 
also assume that “true” image is a blurred version 
of the object we would like to observe, where 
blurring is linear. 7.e. true part of the observed 
image is: 


Prrue(x) = f, Pog (x y)dy (1) 
where g(x,y) is a blurring function that is 
normalised everywhere: 


[gœ y)dy = 1 (2) 


where Ptrye(x) is the true image, po (x) is the 
image we would like to observe. 

Basic image processing and modelling prob- 
lem is that our data have noise and blurred then we 
would like to solve the equation (1) so that as much 
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as possible details in the map can be visualised and 
modelled. This problem can be considered as Wie- 
ner filtering problem (Vega and Rey, 2012). In ge- 
neral, a solution of the equation (1) requires sol- 
ving very large linear systems which is not possible 
in practice. In addition, since we do not know the 
blurring function then the problem becomes even 
more complex; we must find blurring function as 
well as solving the equation (1). When blurring is 
position independent, i.e. g(x,y) = g(x — y) then 
the problem becomes simpler, “true” image is just 
a convolution of the unblurred image with the blur- 
ring function. In this case, it is known that 
Fptur (S) = Fo(s)F (s), i.e. Fourier coefficients of 
the blurred image is the product of Fourier coeffi- 
cients of the unblurred image and that of blurring 
function. 

So, the observation in the Fourier space can be 
written as: 


Fops(S) = Fo(s)F (s) + Fnoise(S) 


Recall that our observations are within the 
nominal resolution sphere |s| < Snom. We assume 
that we know the blurring function - Fj(s), the 
variance and covariance of the noise in the Fourier 
space - Fnoisee We want to know what are the 
details in the image that can be seen (distinguish 
from each other) for a given nominal resolution, 
blurring function and noise level. We will 
approach to this problem step by step. First, we will 
assume that there is no noise Fryojse(s) = 0, no 
blur F} = 1 and the only shortcoming of our obser- 
vations is that we observe only within a resolution 
sphere |s| < Snom, the only problem is the series 
termination effect. After analysis of this problem 
we will move on to assume that there is a position 
independent blurring and we know blurring functi- 
on. We also will assume that the blurring function 
is Gaussian and we know the blurring parameter. 
Then we will consider the cases when there is noise 
and we know the variance and covariance of the 
noise. Then we assume that we do know the shape 
of blurring function but we do not know its para- 
meters (e.g. blurring is Gaussian but we do not 
know variance of the Gaussian). Our problem then 
will be a) what level of details can be seen; b) what 
are the optimal parameters of blurring function. 
Then we will move on and assume that blurring is 
position dependent. However, in this work we will 
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only consider only limited resolution data (series 
termination effect) and position independent Gaus- 
sian blurring. Remaining problems will be dealt 
with in future works. 

To define the resolvability in the scattering 
data we will assume that there are two point sour- 
ces, however because of observational deficiencies 
or intrinsic properties, e.g. mobility of the objects 
observed sources are broadened. We will consider 
two reasons for peak broadening: a) our observati- 
ons are within a certain resolution limit - |s| < 
Snom; b) sources oscillate during observations and 
we observe blurred version of these sources. 

We will define resolvability as follows: With 
a given nominal resolution limit, blurring and noise 
level resolvability is the minimum distance betwe- 
en like point sources that can be distinguished in 
the map. 

We assume that we have two point sources in 
n-dimensional space and they are blurred by a Ga- 


. : . B 
ussian blur with the variance u = m Then the ob- 


served peak for one point source vvould have the 


form: 
Blsl” 


fat) e İL eF eds (a) 
Isl<Snom 

where integral is in n-dimensional space. xs is the 

bilinear form formed with the vectors in image (x) 

and Fourier space (s), i = V—1 is the imaginary 

number. We would like to find minimal non-zero 

root of this equation: 


f:(0) + f(x) — 2f: (=) = 0 (4) 
i.e. we want to find such a distance between point 
sources that image values at the positions of the 
peaks are bigger than in the middle of the vector 
connecting them. 
In the following sections, we will consider 
this equation for one, two and three dimensions 
with and without blurring. 


2. Peak widening and resolvability in 1D 

Let us assume that we have a point source: 
ö(x) in one dimensional case. Its Fourier transfor- 
mation for all s F(s) = 1. Now let us assume that 
we cut the data s € [—Syom, Snom]. Let us calculate 
the inverse Fourier transformation which is the 
following integral: 


“Snom 


= 2Snomjo(2TXSnom) (5) 


2TXSnom 


Therefore, in this case, the above stated 
problem is reduced to the solution of the equation: 


2Snom + 2SnomJo (21x Som) .. 
—4Snomjo(2TXSnom/2) = 0 (6) 
If we let y = Z27TXShnom then the solution of (6) is 
equivalent to finding of the smallest non-zero root 
of the equation 




















1+ jo) 2/00//2) = 0 (7) 
Solution of this equation is Yọ € 4.28 and 
Yo 0.68 
therefore xo = 2755 
4 2 o 2 4 


distance from the centre 


Fig. 1. The graphs of f(x) = jo (2XSnom) function for 
different values of nominal resolution - Snom. This 
graph shows that peaks are broadened as Snom 
decreases. Black line corresponds to Snom = 0.75, blue 
line corresponds to Shom =0.5, red line corresponds to 
Snom = 0.4 and cyan line corresponds to Snom = 0.25. 


Figure 1 demonstrates the effect of nominal 
resolution on peak broadening. As it is expected as 
the nominal resolution increases the calculated 
peaks become sharper. 

If we have two point-sources and observations 
inside the interval |s| < Snom what is the minimal 
distance between peaks so that they can be seen as 
a separate peak. Let us assume that distance 
between point sources is z. Then we can plot the 
following function: 


İo(ZTTXSnom) + jonx S Z)Snom) 
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Fig. 2. The graph of tvvo point-sources for one-dimensio- 
nal case — the function jọ(2TXSnom) + /fo(Zr(x — 
Z)Snom). Resolution is Spom=0.5. Top left: distance 
between sources is 2A, top right: distance between source 
is 1.6A, bottom left: distance between sources 1.36A, 
bottom right: distance between sources is 1.2A 


Fig. 2. demonstrate that inside |s| < 0.5A71 
interval when the distance between peaks are 
decreasing then resolvability power of the map 
become is reduced. When the distance is 2A then 
peaks are clearly visible, when the distance 
becomes 1.6A peaks become less separable, when 
the distance between peaks is 1.36A then peaks are 
marginally separable and when the distance is 1.2A 
then two peaks are merged together and they are no 
longer be seen as two different peaks. It should be 
noted that when the distance is equal to the 
resolvability — 1.36A then, although peaks can still 
be separated they come close to each other - the 
maxima of peaks are no longer at the original peak 
positions. 

Gaussian blurred source with the variance 


B 
u = — has the form: 
8m 


x2 
p(x) =2 ke” (8) 
Fourier transformation of this source is: 
F(s) = e”Bisi”/4 (9) 
Novv vve again assume that the data have been 
observed vvithin the interval se 


İ—Snom? Snom]. Inverse Fourier transformation of 
the Gaussian blurred point source is: 
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fg(x) = — [Pe —B\|s|7/4 ePTXS ds (10) 





Snom 
resulting in 
Vi - (2nx)? — i2mx 
f(x) -2-—e” B Re(erf(—“* - —)) 


VB VB 


VVe can evaluate this ... if vve have 
access to the error function of complex variables. 
This function can be evaluated efficiently using 
Faddeeva function (Poppe and Wijers, 1990) 
defined for complex arguments as: 


w(z) = e~” erfc(—iz) = e7” (1 — erf (—iz)) 


where z is a complex variable and the values of 
this function, in general, are complex numbers. 

Using the relationships between error and 
Faddeeva functions we can write: 





fax) 
Va _(2nx)? 
= 2—_e B 
vB 
vr nom R 27x 
—2-—€ ” 4 Re(e@™*Snomy $ — 
VB vE 
+ (NBsnom ) 


Let us compare the effects of different B values on 
point source broadening. 
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Fig. 3. Point source with blurs. All plots correspond to 
Snom = 0.5, Black line corresponds to B = 0.001A2, 
blue line corresponds to B = 10A”, red line corresponds 
to B=20Å and cyan line corresponds to B=50A”. 


It is clear that as B value increases the peaks 
are broadening thus reducing resolvability of these 
peaks. To observe highly oscillating peaks one 
needs to produce maps with as little as possible 
noise. 
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Fig. 4 shows that for the same resolution and 
distance between Gaussian sources as blurring 
parameter increases peaks become less and less 
resolvable. It should be noted that if there would be 
no noise and blurring parameter would be known 
then we could remove the effect of blurring 
without difficulty by simply using deconvolution. 
The problem becomes more complicated when the 
noise comes into play. All observations in reality 
are with noise, the noise level can be reduced it is 
impossible to remove it completely. 
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Fig. 4. Sum of two Gaussian blurred sources at Snom = 
0.5, dist = 2A for different blurring parameters. Top 
left: B=0.01A~, top right: B-10A7, bottom left B-50A7, 
bottom right B=100A”. 


3. Peak widening and resolvability in 2D 
Now let us assume that we have a point source 
in two dimensions. Its Fourier transformation for 
all s is constant. Now let us assume that we cut the 
data within |s| < Snom. The observed peak will 
have the form: 
fo(x) m ra elds m 
J1(27|x|Snom) 
21|x|Snom (11) 
where J; is the 1* order Bessel function of first 
kind. All other terms are the same as defined 
ht) 
t 
TSzom- If vve uset = 27|x|Syom then the equation 
(4) is reduced to: 
AL) h(/2) 
The solution of this equation is tọ s” 4.7715 
and we obtain that: tọ = 27S, 9, and therefore in 


= 2182 of 


, 1 . 
above. Since lim = 5 we can write fo(0) = 
> 


the absence of blurring resolvability in two- 
2 = 0.7594 —. 


Since these functions have qualitatively the same 
form as in one-dimensional case we do not show 
figures for unblurred cases in two dimensions. 

The observed peak with the blurring parame- 
ter B and nominal resolution Snom is calculated 
using the following integral: 


f(x) = 


dimensional case is: Xo = 








2m Snom 


=f | e~Flsl*/4|5| cos(2m|x||s| cos())ddds 
0 0 


To our knowledge there is no closed form 
expression for this integral, therefore this integral 
is calculated numerically using the statistics 
package R (R core team, 2018) that uses 
QUADPACK for numerical integration (Piessens 
et al., 1983). 


a | 
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Fig. 5. Point source with blur in 2D, Spom = 0.5A7t. 
Black line corresponds to B = 0.001, for blue line B = 
10, for red line B=20 and for cyan line B=50. 


The figure 5 illustrates that for a given 
resolution the observed peaks are broadened as 
blurring increases. 

Again, as in one-dimensional case, let us 
consider the blurred image with two point sources 
and with the observations within the circle |s| < 
Snom - We assume that the distance between these 
point sources is z. Now we can plot the following 
function: 

İBÜZTXSnom) + fg (20(x — Z)Snom) 
to see how peaks widen and become unsolvable as 
blurring parameter increases (figure 6). 
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Fig. 6. The graph of two point sources in two dimen- 
sions with resolution Snom = 0.5A71 and distance bet- 


ween sources d = 2.0A for different blurring parame- 
ters. Top left: B-0.01A7, top right: B=10A7, bottom 
left B-50A 7, bottom right B=100A~. 


The Figure 6 illustrates again that in two di- 
mensions also as blurring parameter — B increases 
peaks become less and less resolvable. 


4. Peak widening and resolvability in 3D 

Now let us consider three-dimensional case 
with point sources. Peaks for which Fourier 
coefficients observed within a nominal resolution 
will have the form: 

foto = 1. eet ds = 
—21|x|Snom COS(27|X|Snom)+sin(2m|x|Snom) _ 

(2n|x\)3 R 
“"iStomİ(Zrilxİsnom) 
2TTShomX 

where j; is the spherical Bessel function of first 
order. All other notations are the same as defined 
above. 


In this case, the solution of equation (4) gives: 


0.8322 : : ; 
Xres = ~ . It is the maximum attainable 
nom 


resolvability in three dimensions in the absence of 
blurring. We will consider blurring with various 
parameters, although in this case the results will be 
qualitatively similar to one- and two-dimensional 
cases. Indeed, it can be seen in figure 7 and 8. 

We next consider the case when the density is 
blurred with a Gaussian blurring function with the 


At 
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parameter B. To add the effect of blurring we again 


need to calculate the integral: 
Bis? 
0 
Isl<Snom 
2 (27|x\|)2 Bs2 . 
= (=) 2 ü = m 4 sin 21|x|Snom 
|x|B 
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Fig. 7. Peak broadening in the case of no blurring 
corresponding to different nominal resolution. These graph 
shows peaks are broadened as Snom decreases. Black line 
corresponds to Snom = 0.75, blue line Snom “0.5, reac 
line Shom = 0.4 and cyan line Shom = 0.25. 








Peak Height 
Peak Height 












































2 0 2 4 2 0 2 4 
Distance Distance 
S 
oo x 
əni s 
E o E 
ə 37 $ 34 
I I 
x + | xo 
$ s $ so] 
a a 
a = 
~j s 1 
24 e| 
m T T T T S T T T T 
2 0 2 4 2 0 2 4 
Distance Distance 


Fig. 8. Sum of two peaks for three-dimensional case 
with Gaussian blurring. Top left corresponds to 
B-0.001A 7, top right corresponds to B=10A~, bottom 
left corresponds to B =50Å? and bottom right 
corresponds to B=100A* 
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All terms and variables are as defined before. 
The first term on the right side of this equation is 
the peak if there is no series termination and the 
second terms accounts for the series termination. 

Again, we need to plug in this function into 
the equation (4) and solve it. As expected as B 
values increase the solution of this equation gives 
larger and larger distances meaning that when 
amplitude of blurring increases the resolvability of 
peaks decreases. 

Again, the figure 8 shows that as blurring 
amplitude increases peak separability and therefore 
resolvability decreases. Detailed analyses of cases 
with noise and anisotropic mobility will be the 
subject of future study. 


Conclusion and Future Perspectives 

In this work, we presented an approach to the 
problem of resolution and resolvability for scatte- 
ring methods used in macromolecular structural 
biology. We focused on the cases without noise. 
We assumed that there are only two shortcomings 
of the observations: 1) diffraction limit, i.e. obser- 
vations are made within the nominal resolution 
İsİ < Snom; 2) sources for which we are doing ob- 
servations are oscillating and as a result we observe 
a blurred version of the object. We considered both 
cases for one-, two- and three-dimensional cases. 
We showed that as resolution decreases peaks get 
broader and thus becoming less resolvable as expec- 
ted. We also found a closed form relationship bet- 
ween the nominal resolution and peak resolvability 
for one-, two- and three-dimensional cases. Similar 
effects are observed when objects are blurred with a 
Gaussian function — peaks get broader as the ampli- 
tude of oscillation increases and they become less 
resolvable. It has already an implication to structural 
biology that signal level in different parts of the ob- 
served density will be different and will correspond 
to the amplitude of oscillation of the molecule cor- 
responding to these parts of the density. 

Future work will focus on a theoretical treat- 
ment of the effect of noise on the peak resolvabi- 
lity. Ultimate purpose of this work, that will be de- 
alt in future works, is to derive an algorithm to cal- 
culate position and orientation dependent resolva- 
bility in the density corresponding to noisy data 
with anisotropic signal strengths. The planned 
work will answer to one of the long-standing ques- 
tions in structural biology: what are the details that 


can be observed in the map if we do know the no- 
minal resolution, the noise level, relative and abso- 
lute amplitudes of mobility for different parts of 
the calculated (observed) maps. The result of this 
work will help the unification of resolvability used 
in different scattering methods. However, we wo- 
uld like to stress that our main concern is macro- 
molecular crystallography and single particle cryo- 
electron microscopy. 
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Bir-, iki və üçölçülü fəzalarda ayırdetmə və ayırdetmə qabiliyyəti 
G.M. Qasimova!”, R. Ç. Masmoliyeva’, Q. N. Mürşüdov” 


" AMEA Molekulyar Biologiya və Biotexnologiyalar İnstitutu 
?Azərbaycan Dövlət Neft və Sənayə Universiteti 
"TTM Molekulyar Biologiya Laboratoriyası, Kembric, Böyük Britaniya 


Bu məqalə difraksiya məhdudluğu və bulanıqlıq mövcud olan hallarda səpilmə eksperimental metodlarının 
(məsələn, kristal diffraksiyası, elektron mikroskopiyası) ayırdetmə və ayırdetmə qabiliyyəti probleminə bir 
yanaşmaya həsr olunur. Birölçülü, ikiölçülü və üçölçülü halların hər biri ayrı-ayrılıqda nəzəri analiz edil- 
mişdir. Nominal ayırdetmə bulanıqlıqdan asılı olan zirvə ayırdetmə qabiliyyətinin ifadələri alınmış və təd- 
biq olunmuşdurlar. Biz difraksiya məhdudluğu və bulanıqlıq böyüdükcə nöqtə mənbələrinin bir-birindən 
ayrı görünməsinin çətinləşdiyini göstərdik. 


Açar sözlər: Saflaşdırılma, elektron kriomikroskopiyası, Furye qabığının korrelyasiyası, Furye çevirməsi, 
Qauss paylanması 


Pazpemenne n pa3peHIHMOCTPE B O/IHOM, /IBYX H TpeX H3M€peHHSsX 
T.M. TacbıoBal?, P./Pk, MacmalıneBal, T.H. MypmyaoB!? 


! Hucmumym monekynapnoŭ 6uonoeuu u 6uomexnonoeutt HAH Asepöaüğərana 
? Asepöaüərcanckxulü socyöapcmeennəltü nepmanoŭ u npoxpiulnenHətü yuusepcumem 
İMRC /Ta6öopamopus monrekynapnoŭ Guonozuu, Kem6pudoic, Benuxo6pumanusı 


B ərofi cTaTbe OMMCbIBacTCA MOZXOA K MmpoOJIeMe pa3pelieHHA H pa3peHHMOCTH B M€TO/TaAX paccesHHii 
(HarıpHMep, MAPpaKUHA peHTTEHOBCKHX /IyHCİİ, 3IeKTPOHHad MUKPOCKOMNMA) TIpH HANHYHH paspbiBa pana 
@Mypbe H pa3sMbITHA. ÖHHO-, HByX- H TpEXMEPHBI€ c/yuaH paccMaTpuBatoTca oTZeubHO. İloyryueHbi H 
rpoaHasır3HpOBaHBI (ÞOpMyJIbI, CBA3bIBAIOINHE ƏİPþekKTHI HOMHHAIbHOTO pa3peleHHA H pa3MbITHA C 
MaKCHMAJbHOÑ pa3pelarorei CHOCOÕHOCTÞIO. MbI NOKa3bIBAEM, YTO KaK pa3MBITHE, TAK H pa3pbIB pana 
Dypbe paciupsAoT IHKA TOHCHHOTO HCTOYHHKA, TEM CAMbIM CHWKAIOT HX paspernaroIyIo CT1OCOÖHOCTB. 


Kunoueevpie croea: Ymounenue, 31ekmponnaa Kpuo-mukpockonusa, Kopperayua o6onouek Pypve, npe- 
oöpasosanue Pypve, pacnpedenenue laycca 


36 


